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PATENT 

Attorney Docket No.: 091 03-0 16220US 
PULSE OXIMETER SENSOR WITH PIECE-WISE FUNCTION 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application is a continuation of US Application No. 09/836,050, filed 

5 April 16, 2001, now US Patent No. , which claims the benefit of US Provisional 

Application No. 60/198,109, filed April 17, 2000, which disclosures are incorporated 
herein by reference. 

BACKGROUND OF THE INVENTION 
10 The present invention relates to oximeter sensors having a memory. 

Pulse oximetry is typically used to measure various blood flow 
characteristics including, but not limited to, the blood-oxygen saturation of hemoglobin in 
arterial blood, and the rate of blood pulsations corresponding to a heart rate of a patient. 
Measurement of these characteristics has been accomplished by use of a non-invasive 
15 sensor which passes light through a portion of the patient's tissue where blood perfuses 
the tissue, and photoelectrically senses the absorption of light in such tissue. The amount 
of light absorbed is then used to calculate the amount of blood constituent being 
measured. 

The light passed through the tissue is selected to be of one or more 
20 wavelengths that are absorbed by the blood in an amount representative of the amount of 
the blood constituent present in the blood. The amount of transmitted or reflected light 
passed through the tissue will vary in accordance with the changing amount of blood 
constituent in the tissue and the related light absorption. For measuring blood oxygen 
level, such sensors have been provided with light sources and photodetectors that are 
25 adapted to operate at two different wavelengths, in accordance with known techniques for 
measuring blood oxygen saturation. 

Various methods have been proposed in the past for coding information in 
sensors, including pulse oximeter sensors, to convey useful information to a monitor. For 
example, an encoding mechanism is shown in Nellcor U.S. Patent No. 4,700,708. This 
30 mechanism relates to an optical oximeter probe which uses a pair of light emitting diodes 
(LEDs) to direct light through blood-perfused tissue, with a detector picking up light 
which has not been absorbed by the tissue. The operation depends upon knowing the 
wavelength of the LEDs. Since the wavelength of LEDs can vary from device-to-device, 
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a coding resistor is placed in the sensor with the value of the resistor corresponding to the 
actual wavelength of at least one of the LEDs. When the oximeter instrument is turned 
on, it first determines the value of the resistor and thus appropriate saturation calculation 
coefficients for the value of the wavelengths of the LEDs in the probe. 
5 Other coding mechanisms have also been proposed in U.S. Patent Nos. 

5,259,381; 4,942,877; 4,446,715; 3,790,910; 4,303,984; 4,621,643; 5,246,003; 3,720,177; 
4,684,245; 5,645,059; 5,058,588; 4,858,615; and 4,942,877, the disclosures of which are 
all hereby incorporated by reference. The '877 patent in particular discloses storing a 
variety of data in a pulse oximetry sensor memory, including coefficients for a saturation 

10 equation for oximetry. 

Nellcor pulse oximeter sensors are encoded with a resistor (RCAL) value 
that corresponds to the wavelength(s) of the LED(s) within the emitter, such as described 
in Pat. No. 4,700,708. Nellcor pulse oximeter instruments read this resistor coding value 
and use it as a pointer to a look-up table that holds the proper set of coefficients for that 

15 sensor for calculating arterial oxygen saturation (SpOj). The function that converts the 
measured red and IR signal modulation ratio R (also known as the "ratio of ratios" or 
"rat-rat") to a calculated saturation value is derived from the basic form of the Lambert- 
Beer Law: 

R _ \n(I ] /I 2 ) red _ S-/3g„ b +(l-S)-/3% _ Sc l+ (l-S)-c 2 
ln(/,// 2 ),> S-/3Z 2Hb +(l-S)-j3Z b S-c 3 +(l-S)-c 4 

20 where and L refer to detected light signals at two different points in the cardiac cycle, 
and the /?'s refer to the characteristic light absorption properties of oxygenated and 
deoxygenated hemoglobin. When solved for the saturation (5), the result takes on the 
form: 

S P 0 2 =5-100- c 2 -c 4 R 10Q 

(c 3 -c 4 )tf + (c 2 -c x ) 

25 (2) 
Equation 2 can be further simplified to require only three constants (by, for 
example, dividing each constant by c 2 ), but will be used as shown for the remainder of 
this description. Although theoretically based, the four constants C1-C4 are empirically 
determined. Theoretical values for the constants are insufficient primarily due to the 

30 complexities of light scattering and sensor optics. The values of the sets of constants (c/ 
through c 4 ) vary with each resistor coding bin (each "bin" corresponding to a range of 
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different characterized LED wavelengths). Multiple sets of coefficients (bins) are 
provided within a lookup table in Nellcor oximeters. When calculated Sp0 2 values 
according to Eq. 2 are less than 70%, a revised value of Sp02 using a linear function is 
used: 

5 Sp0 2 =c 5 -c 6 R, (3) 

where both c$ and c$ vary with the resistor coding value. This linear function was found 
to better match SpC>2 (arterial oxygen saturation as measured by a pulse oximeter) with 
Sa02 (the true value of arterial oxygen saturation, as measured directly on a blood 
sample) in observations made at low saturations. 
10 A limitation of this method is that the proper calibration of the pulse 

oximetry sensor can be accomplished only if the relationship between the signal 
modulation ratio (R) to blood Sa0 2 conforms to one of the pre-encoded sets of calibration 
coefficients. 

A further limitation of this method is that the relationship between R and 
15 Sa0 2 of the pulse oximetry sensor may not be linear in a low-saturation region, or that the 
breakpoint may not optimally be located at 70% Sp0 2 . 

A yet further limitation of this prior art method is that the functional 
relationship between the true arterial oxygen saturation and the measured signals may not 
fit a single function over the entire span of the measurement range. 

20 

SUMMARY OF THE INVENTION 
The present invention takes advantage of a memory in the sensor to 
provide enhanced performance. In one embodiment, not only are the sensor's specific 
calibration coefficients stored in a memory in the sensor for the formula to determine 
25 oxygen saturation, but multiple sets of coefficients are stored. The multiple sets apply to 
different ranges of saturation values to provide a better fit to occur by breaking the R to 
Sp02 relationship up into different pieces, each described by a different function. The 
different functions can also be according to different formulas for determining oxygen 
saturation. 

30 In another aspect of the invention, the sensor can store a variable 

breakpoint between the two functions used for oxygen saturation. The two functions 
could either be separate formulas or the same formula with different coefficients. This 
allows optimization to a value other than the 70% breakpoint of the prior art. 
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In another aspect of the present invention, the sensor can store more than 
one breakpoint to create more than two functions describing the R to Sp02 relationship. 

In yet another aspect of the present invention, a spline function is used, 
breaking up the R to Sp02 relationship into an arbitrary number of regions. 
5 In one embodiment, the coefficients stored in the sensor memory 

correspond to a non-linear curve for low saturation values below 70% or some other 
breakpoint(s). 

Each of the methods described here improve the fit between the chosen 
mathematical function and the arterial oxygen saturation by breaking the relationship into 
10 subsets of the full measured range and determining optimum coefficients for each range. 
Spline-fitting, in this context, similarly breaks the full measurement range into subsets to 
efficiently describe the numerical relational between the underlying tissue parameter of 
interest and the actual signals being used to estimate its value. 

For a further understanding of the nature and advantages of the invention, 
1 5 reference should be made to the following description taken in conjunction with the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a block diagram of a pulse oximeter system incorporating the 
20 present invention. 

Fig. 2 is a graph of R (signal modulation ratio) versus oxygen saturation 

(Sa0 2 ). 

Fig. 3 is a diagram of the contents of a sensor memory according to the 

invention. 

25 Fig. 4 is a graph of oxygen saturation versus R to illustrate the 

embodiment for spline or curve fitting to a predefined set of knots. 

Figures 5A, 5B, 6A and 6B are graphs illustrating the improved curve 
fitting of the embodiments of the invention versus the prior art. 

30 DESCRIPTION OF THE SPECIFIC EMBODIMENTS 

Sensor Reader/Monitor 

Fig. 1 is a block diagram of one embodiment of the invention. Fig. 1 
shows a pulse oximeter 17 (or sensor reader) which is connected to a non-invasive sensor 
15 attached to patient tissue 18. Light from sensor LEDs 14 passes into the patient tissue 
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18, and after being transmitted through or reflected from tissue 18, the light is received by 
photosensor 16. Either two or more LEDs can be used depending upon the embodiment 
of the present invention. Photosensor 16 converts the received energy into an electrical 
signal, which is then fed to input amplifier 20. 
5 Light sources other than LEDs can be used. For example, lasers could be 

used, or a white light source could be used with appropriate wavelength filters either at 
the transmitting or receiving ends. 

Time Processing Unit (TPU) 48 sends control signals to the LED drive 32, 
to activate the LEDs, typically in alternation. Again, depending on the embodiment, the 

1 0 drive may control two or any additional desired number of LEDs. 

The signal received from input amplifier 20 is passed through two 
different channels as shown in the embodiment of Fig. 1 for two different wavelengths. 
Alternately, three channels for three wavelengths could be used, or N channels for N 
wavelengths. Each channel includes an analog switch 40, a low pass filter 42, and an 

1 5 analog to digital (A/D) converter 38. Control lines from TPU 48 select the appropriate 
channel at the time the corresponding LED 14 is being driven, in synchronization. A 
queued serial module (QSM) 46 receives the digital data from each of the channels via 
data lines 79. CPU 50 transfers the data from QSM 46 into RAM 52 as QSM 46 
periodically fills up. In one embodiment, QSM 46, TPU 48, CPU 50 and RAM 52 are 

20 part of one integrated circuit, such as a microcontroller. 

Sensor Memory 

Sensor 15, which includes photodetector 16 and LEDs 14, has a sensor 
memory 12 associated with it. Memory 12 is connected to CPU 50 in the sensor reader or 

25 monitor 17. The memory 12 could be packaged in a body of the sensor 15 or in an 

electrical plug connected to the sensor. Alternatively, the memory 12 could be packaged 
in a housing which is attachable to an external surface of the monitor or the memory 12 
could be located anywhere in a signal path between the sensor body and the monitor. 
Specifically, according to some preferred embodiments, a content of the sensor memory 

30 12 could be constant for all sensors associated with a particular sensor model. In this 
case, instead of putting an individual memory 12 on each sensor associated with this 
model, the memory 12 could instead be included in a reusable extension cable associated 
with the sensor model. If the sensor model is a disposable sensor, in this case a single 
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memory 1 2 would be incorporated into a reusable extension cable. The reusable cable 
could then be used with multiple disposable sensors. 

Fig 2 is an example of a graph of the ratio of ratios (R) on the X axis 
versus oxygen saturation (SaC>2) on the Y axis. Shown is a breakpoint 52. In the prior 
5 art, a breakpoint of 70% was pre-defined in the monitor software. To the right of the 
breakpoint (oxygen saturations between 70-100%) a formula was used with four 
coefficients. To the left of the breakpoint in the prior art, a linear equation was used with 
two coefficients. The present invention provides increased flexibility and accuracy by 
using a non-linear formula for the portion of the curve to the left of breakpoint 52. By 

10 using a memory chip in the sensor itself, it is possible to actually store these coefficients 
on the memory chip, as well as the separate coefficients for the higher saturation values. 

In another embodiment of the invention, breakpoint 52 can be stored in the 
memory chip, and chosen to optimize the curve fitting for the two sets of coefficients. In 
other words, a better fit to the two curves may be obtained if the breakpoint is 68%, for 

15 example. In an alternate embodiment, multiple breakpoints and curves might be used. In 
addition, rather than using the same formula, different formulas could be used for 
different sections in another embodiment. 

Fig. 3 illustrates the contents of sensor memory 12 of Fig. 1. As 
illustrated, in a first section of memory 54 are stored a first set of coefficients. A second 

20 portion of memory 56 stores a second set of coefficients. Finally, in a third section of 
memory 58, the breakpoint 52 is stored. Different combinations of these elements could 
be stored in different memories. For example, the breakpoint could be left out of some, 
and in others a breakpoint may be provided with only one set of coefficients (with the 
other set of coefficients in the monitor). Alternately, a breakpoint might be implied from 

25 a sensor model number which is stored in the memory, or some other ED value. 



yg-equation : 

In one embodiment, an enhanced form of the curvilinear function is used. 
Instead of using Eq.3 (linear) in the lower saturation region, Eq.2 (non-linear) is used for 
30 both the upper and lower saturation regions. The breakpoint that defines when to switch 
coefficients from an upper-region set to a lower-region set is defined by another 
coefficient. The breakpoint can be programmed either as a value of i?, or as a value of 
SpC>2. With the breakpoint defined as a value of/?, the algorithm becomes: 
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(c-d)R + (b-a) 



100- 



\R<c s : a = c { ,b — c 2 ,c - c 3 ,d = c 4 
I R > c s :a = c 6 ,b = c n ,c = c z ,d = c 9 



(4) 



Curve Fitting 

Curve fitting to multiple regions follows the same methodology as fitting 
5 to a single region. Simply put, the data is partitioned into separate regions and 

coefficients are determined for each region separately. Commercially available software 
programs are available, (for example, Mathcad, (Mathsoft, Inc., Cambridge, MA). The 
process can also be found in, for example, Data Reduction and Error Analysis for the 
Physical Sciences (Philip Beviyton, McGraw-Hill, New York 1969, Ch. 1 1 - Least 
1 0 squares fit to an arbitrary function). 



Spline Fitting 

An alternate embodiment uses either spline (curve) fitting, or linear or 
higher order interpolation to a predefined set of Sp0 2 vs R values ("knots"). A "knot" is a 
1 5 term of art in spline fitting that refers to an x-y pair corresponding to a node on a line, 
with a number of such knots defining the line. Spline fitting is a technique for 
interpolation. 

For instance, the values oiR at specifically defined SpC>2 values would be 
stored in the sensor memory. An example of this looks like: 
20 R = a b c 

SpO 2 = 100 95 90 

Alternatively, though less preferably, the independent variable could be swapped: 

R- 0.5 0.6 0.7 

Sp0 2 = x y z 
25 a) Only the bold values (e.g., a, b and c) would need to be stored with fixed, pre- 
selected spaced values of Sp0 2 (equally spaced or unequally spaced). Or, 
alternatively, preselected values of R. 

b) An alternative approach would store within the sensor memory the 
iS/7(92(minimum) and Sj902(maximum) values of the spline range, the number of 

30 knots that will be defined, and the sequence of defined values of R for those knots. 

c) A further alternative approach could store both SpC>2 and the associated R value 
for each knot. 
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For each of these options, the instrument would use a spline-fitting algorithm, preferably 
a cubic spline, to determine the Sp02 at the measured value of R according to the stored 
values (an alternative could be a linear or higher order interpolation algorithm). 

Fig. 4 illustrates the cubic spline method. Fig. 4 is a graph of oxygen 
5 saturation vs. R for a particular sensor emitter. Thus, instead of storing the coefficients as 
in the prior art method, the actual R or oxygen saturation values are calculated and stored 
in the sensor memory for that particular sensor's characteristics (e.g., emitter 
wavelengths). When the oximeter measures the signal level of the light detector, it 
determines an oxygen saturation value by determining the point on the curve associated 

10 with the calculated R value between two of the sample points shown in Fig. 4. 

There exists a trade-off in the number of knots defined and the amount of 
memory required to store them. Too few knots requires very little storage memory, but 
may not adequately describe the functional relationship; too many over-defines the curve 
and consumes more memory. The inventors have found that knots spaced 5%-10% apart 

1 5 give adequate results. 

Cubic Spline Calculation : 

The process for cubic spline interpolation is known to those skilled in the 
art. Intrinsic in using the spline method is that the value of R needs to be determined first 

20 before being translated to SpC>2. The preferred process for spline interpolation can be 
accomplished using the functions provided in Mathcad, and treats the endpoints with 
cubic functions. Other references for cubic spline interpolations are available. 

The process of finding the coordinates of the knots in empirical data with a 
significant amount of noise may require an additional step. Commercially available basic 

25 curve fitting programs may be used (sigmaPlot, or TableCurve, or Mathematical for 

instance) to determine a best-fit functional approximation to the data. Alternately, one 
can perform a least-squares fit of an arbitrarily chosen analytical function and pick the 
values of R at the knot locations (SaOi values). The analytical function can be an 
overlapping piece-wise polynomial (e.g., linear or parabolic), or the curvilinear equation 

30 of Eq. 1 or Eq. 4. Another approach is to perform a least-squares selection of the knots 
directly. 

Figure 5A shows the conventional curve fitting of the prior art, wherein a 
linear relationship is used below 70% saturation, with a curvilinear approach above 70%. 
The residual error due to an imperfect fit to the actual R to Sa0 2 response for the 



curvilinear approach above 70% saturation is illustrated by curve 60, while the residual 
error of the linear interpolation approach below 70% is illustrated by dots 62. Figure 5B 
illustrates the use of curvilinear fits in both regions, with a different curvilinear curve 64 
being used below 70%. In this instance, a much improved fit is provided. In both 
5 figures, the smaller dotted line 66 corresponds to the use of a single curvilinear fit across 
both regions, which is also not as accurate, having a much higher error characteristic 
compared to the curves of the invention, 64 and 60 of Figure 5B. 

Figure 6 A and 6B show a plurality of knots as circles 70 on the graphs. 
Dotted line 72 of Figure 6A illustrates a linear interpolation fit to these knots, which 

10 shows a residual error prone result with multiple loops. In Figure 6B, on the other hand, 
the present invention using a cubic spline fitting approach provides a dotted line 74 which 
is a more accurate fit to the knots 70. 

As will be understood by those of skill in the art, the present invention 
may be embodied in other specific embodiments without departing from the essential 

15 characteristics thereof. For example, any function can be used for the formulas for 

determining oxygen saturation, not just the ones described. For a limited sensor memory, 
the function representation may be compressed. Any representation of a function could 
be used. Calibration coefficients may be based on more or different characteristics than 
the sensor's LED wavelength(s). For example, other LED emitter characteristics or 

20 sensor design characteristics can be factors in the sensor's calibration coefficients. 

Additionally, the formula for calculating oxygen saturation may be a 
function of more than the ratio of ratios; for example, other input variables such as signal 
strength, light levels, and signals from multiple detectors could be used. 

This methodology for piece-wise fitting is not limited to oximetry. This 

25 method is useful when the relationship between the measured signal and reference value 
observed during calibration is not adequately described by a single function or set of 
coefficients over the whole measurement range. The relationship may be broken into 
subsets, and a piece-wise continuous set of functions may be used to describe the 
relationship. For example, other blood or tissue constituents could be calculated, such as 

30 carboxyhemoglobin, methemoglobin, bilirubin, glucose, lactate, etc. Accordingly, the 
foregoing description is intended to be illustrative, but not limiting, of the scope of the 
invention which is set forth in the following claims. 
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